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Temperature fluctuations in an atmospheric convective boundary 
layer are investigated by means of Large Eddy Simulations (LES). 
A novel statistical characterization for both weak temperature fluc- 
tuations and strong temperature fluctuations has been found. De- 
spite the nontriviality of the dynamics of temperature fluctuations, 
our data support the idea that the most relevant statistical properties 
can be captured solely in terms of two scaling exponents, character- 
izing the entire mixed layer. Such exponents control asymptotic (i.e. 
core and tails) rescaling properties of the probability density functions 
of equal-time temperature differences, A r 6, between points separated 
by a distance r. A link between statistical properties of large tem- 
perature fluctuations and geometrical properties of the set hosting 
such fluctuations is also provided. Finally, a possible application of 
our new findings to the problem of subgrid-scale parameterizations for 
the temperature field in a convective boundary layer is discussed. 
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Abstract 



1 Introduction 



Temperature in an atmospheric boundary layer (ABL) is typically convected 
by a velocity background, v, and diffuses by the action of molecular motion 



and/or small-scale turbulent eddies. The basic equation governing such a 
process is the well-known advection-diffusion equation (Pielke, 1984) for the 
(potential) temperature, 9, 

d t 9 + v a d a 6 = D d 2 9 + S e , (1) 

where Sg represents the sources and sinks of heat, eventually present within 
the domain, v is the velocity field advecting the temperature, and D may 
represent either the diffusion coefficients or, alternatively, an eddy diffusion 
coefficient if one intends to focus on the large scale behavior (i.e. large eddies) 
of 9 and thus needs to parameterize, in some way, small scale temperature 
dynamics. Repeated indexes are summed. 

Here, we use the following short notations: d t = d • /dt; 9* = d • /dxi, 
i = 1, ■ • • ,3; d 2 = E?=i 

In many situations, temperature dynamics, driven by the velocity field via 
the advection term in Eq. (0), does not react back on the velocity field. This 
is the case, for instance, in a neutrally stratified boundary layer (occurring, 
e.g., under windy conditions with a complete cloud cover) where buoyancy 
forces are negligible compared to the other terms in the Navier-Stokes (NS) 
equations ruling the velocity field dynamics. In this case, temperature be- 
haves as a passive scalar as a good approximation. 

As a matter of fact, during the diurnal cycle neutral stratification is rarely 
observed in the ABL (Garratt, 1999). More frequently, one observes stable 
conditions (occurring, e.g., at night in response to surface cooling by long- 
wave emission to space) or unstable, convective, ones (occurring, e.g., when 
strong surface heating produces thermal instability or convection). The role 
of temperature is active in both cases, that amounts to say that temperature 
drives the velocity field dynamics through the buoyancy contribution (usu- 
ally modeled by means of the well-known Boussinesq coupling). The latter 
contribution is now the leading one in the NS equations. 

The main feature characterizing both active and passive scalar dynamics 
is the presence of strong fluctuations in the temperature field. Such fluctua- 
tions affect the whole range of scales involved in the temperature dynamics, 
from the largest scales of motion to the smallest ones where diffusive effects 
become important. This huge excitation of degrees of freedom, gives mean- 
ing to the term "scalar turbulence" recently used to denote the dynamics of 
temperature fluctuations (Shraiman and Siggia, 2000). 

The first unpleasant consequence of persistent fluctuations is the fail- 
ure of any attempt to construct dimensional theories for the statistics of 
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temperature fluctuations (Frisch, 1995), quantitatively defined as tempera- 
ture differences A r 9 between points separated by r: A r 9 = 8(r,t) — 9(0, t). 
The common strategy for dimensional approaches consists to define typical 
length/time scales and typical amplitude for the fluctuations of the unknown 
fields (e.g., A r 9 and A r v) and then try to balance the various terms in the 
basic equations (e.g. Eq. ([]]) coupled to the NS equations) to deduce predic- 
tions for A r 9 and A r v as a function of the separation r. 
This is, for instance, the essence of the first dimensional theory for scalar 
turbulence due to Kolmogorov, Obukhov and Corrsin (1949) and Bolgiano 
(1959). As a result of these theories, probability density functions (pdfs), 
P(A r 9), of temperature differences, A r 9, obey the following simple rescal- 
ing: P(A r 9) = r~ a P(A r 9/r a ), where P is a function of the sole ratio A r 9/r a . 
Such property immediately implies that, dimensionally, one has A r 9 ~ r a 
and, for the p-th moment of A r 9: ((A r 9) p ) ~ r^>, with ( p (known as scaling 
exponents) linear function of p: ( p = ap 

The linearity of ( p vs p reflects the fact that only one parameter, a, is nec- 
essary to explain most of the statistical properties of 9. One has, in other 
words, 'single-scale fluctuations'; that amounts to say that it is irrelevant 
which part of the probability density function of temperature differences is 
sampled to define a typical fluctuation. 

Rather than the above simple scenario predicted by dimensional theories, 
turbulent systems show an infinite hierarchy of 'independent' fluctuations 
(Frisch, 1995), that amounts to say the strong dependence on the order p 
considered to define the typical fluctuation. More quantitatively, turbulent 
systems exhibit a nonlinear behavior of C P vs where the infinite set of expo- 
nents, ( p , select different classes of fluctuations. The departure of the actual 
curve of £ p vs p from the linear (dimensional) prediction is named "intermit- 
tency" or "anomalous scaling" (Frisch, 1995). Intermittency is probably the 
most representative property characterizing a turbulent system. 

Our aim here is to provide a statistical characterization of temperature 
fluctuations in a convective boundary layer dominated by well-organized 
plumes, simulated by a large eddy simulation model (Moeng, 1984). 
We shall focus, in particular, on the statistical characterization of two differ- 
ent classes of fluctuations: weak temperature fluctuations, mainly occurring 
in the inner plume regions and strong temperature fluctuations, associated 

1 Such curve must be concave and not decreasing, as it follows from general inequalities 
of probability theory (see, e.g., Frisch (1995)). 
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to the plume interfaces (see Fig. 1). As we shall show, weak fluctuations are 
associated to linear behavior of scaling exponents ( p vs p for small p's, while 
strong fluctuations (captured by large p's) cause the so-called intermittency 
saturations, i.e. the curve ( p vs p tends a to a constant value, Coo, for p large 
enough. The saturation exponent Coo is simply connected (see Sec. |5.1.1|) 
to the fractal dimension of the set hosting large temperature excursions: 
Coo = d — D F , where d and D F are the usual dimension of the space and the 
fractal dimension of the large temperature fluctuation set, respectively. 
Despite the complexity of temperature fluctuation dynamics in a convective 
boundary layer (CBL), reflected in the strong intermittency of the system, 
only two exponents are necessary to capture most of the statistics of tem- 
perature fluctuations. 

It is worth emphasizing that the same statistical characterization has been re- 
cently found for two dimensional idealized models of both passive (see Frisch 
et al (1999), Celani et al (2000) and Celani et al (2001a)) and active scalar 
turbulence (see Celani et al (2001b)) simulated by means of both direct nu- 
merical simulations and Lagrangian methods (see, e.g., Frisch et al (1998)). 
This points toward the generality of our new findings within the context of 
scalar transport. 

2 Statistical tools 

The aim of this section is to provide a quick summary of the statistical tools 
we have exploited to characterize temperature fluctuations in a CBL. 
The basic and well-known indicator is the probability density function, P(Ar-xd), 
of temperature differences, A r . x 0, over a scale r, defined as: 

A r - x = 9(x + r,t)-e(x,t) . (2) 

The pdf P(Ar-x&) will depend other than on r also on x if the system is not 
homogeneous. In our CBL, we have homogeneity along x-y planes, but not 
along the vertical direction z. We thus shall have a dependence on the vertical 
coordinate z. Moreover, in the analysis of Sec. |5|, separations r will be taken 
along x-y planes in the direction forming an angle of 7r/4 with the geostrophic 
wind direction. We shall thus denote our pdf simply as P(A r - z 6). The choice 
for the direction 7r/4 has been done in order to reduce the contamination 
of the scaling exponents by anisotropic effects. More details on this 'magic' 
angle can be found in Celani et al. (2001a). 
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By definition, weak fluctuations (i.e. small values of \6{x + r,t) — 6(x, t) | 
with respect to a typical fluctuation defined as = [(9 2 ) — (O) 2 } 1 ^ 2 ) are 
associated to the pdf core while, on the contrary, large fluctuations (i.e. 
\9(x + r,t) — 8(x, t)\ 3> cr^) are associated to the pdf tails. 
The above considerations can be easily paraphrased in terms of the moments, 
S p (r;z) = ((A r . z 9) p ), known as structure functions. Large fluctuations are 
captured by large p's, while weak fluctuations by small p's. 

Let us now introduce two possible behaviors for the pdf P(A r . z 6). As we 
shall show in the sequel, such behaviors will characterize, within the entire 
mixed layer, the statistical properties of weak temperature fluctuations and 
strong temperature fluctuations, respectively. 

Self-similar behavior 

In terms of probability density functions of A r - Z 6, such a behavior is 
defined by the rescaling property: 

P(M=^(^) ■ (3) 
It can be immediately verified from the definition of moments: 

/+oo 
p{A r , z e){A r , z eyd{A r , z e) (4) 
-oo 

that (0) is equivalent to the following behavior for the structure functions, 

S p {r; z) ~ A Z) with C { P Z) = a (z) P, ® {z) > , (5) 

that is a linear behavior, with the factor a'^ in general depending on the 
elevation z within the mixed layer. 

It is worth noticing that (^]) and (|5[) does not necessarily imply a Gaus- 
sian shape for P(A r - z 6). On the contrary, if P(A r . z 9) is Gaussian then 
(H) (and thus (||)) are immediately satisfied as it follows from the well- 
known property of Gaussian statistics (see, e.g., Frisch (1995)): S-z p (r; z) = 
(2p — 1)!! [52 (r; z)] p , together with the assumption that 52(r; z) ~ r 2a(z) . 

Intermittency saturation (i.e. the strongest violation of dimensional pre- 
dictions): 
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In terms of P(A r . z 6), intermittency saturation is defined as 



r C«? / a a \ 

P{A r , z 6) = -^Q i -g-J for \A r , z 6\>Xa^ (A > 1) , (6) 

where Q is some function (not determined a priori) which does not depend 
on the separation r. 

In terms of cumulated probabilities, i.e. the sum (integral) of the pdfs over 
the large temperature fluctuations (i.e. for |A r;z 0| > \a^ z \ with A > 1), 
defined as: 

/— Aa-M r+oo 
P(A r , z 9)d(A r . z 9) + / P(A r . z B)d(A r;z 6), 
-oo J Ao-M 

(7) 

saturation is equivalent to the following power law behavior, holding for 
different values of A > 1: 

Prob[\A r . z 9\ > Xa {z) ] ~ r c ~ . (8) 

The scaling exponents, Coo , can be thus easily extracted by measuring the 
slope of log {Prob[\A r]Z 9\ > Acr^]} vs logr. 

Finally, it is worth observing that, in terms of structure functions inter- 
mittency saturations means: 



S p (r;z)~ r & } with C£° = C£\ for p > p m t , (9) 

as one can easily verify from the definitions of moments (Q) and from (|6|). 
The scaling exponents (p^ thus tend to a constant value Coo'* for orders, p's 
large enough. Such behavior justifies the word 'saturation' to denote the laws 
d) and (|). 



3 The Large-Eddy simulation model 

In order to gather statistical informations on the turbulent structure of a 
CBL, we used the LES code described in Moeng (1984) and Sullivan et al 
(1994). Such model has been widely used and tested to investigate funda- 
mental problems in the framework of boundary layers (see, e.g., Moeng and 
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Wyngaard (1989), Moeng et al. (1992), Andren and Moeng (1993), Moeng 
and Sullivan (1994), among the others). 

For this reason we confine ourselves only on general aspects of the LES strat- 
egy Details can be found in the aforementioned references. 

The key point of the LES strategy is that the large scale motion (i.e. 
motion associated to the large turbulent eddies) is explicitly solved while 
the smallest scales (typically in the inertial range of scales) are described 
in a statistical consistent way (i.e. parameterized in terms of the resolved, 
large scale, velocity and temperature fields). This is done by filtering the 
governing equations for velocity and potential temperature by means of a 
filter operator. Applied, e.g., to the potential temperature field 9, the filter 
is defined as the convolution: 

»(x) = J »{x')G{x - x')dx' ,10) 

where 9 is the filtered variable and G(x) is a tridimensional filter function. 
The field 9 can be thus decomposed as 

9 = 9 + 9'. (11) 

Applying the filter operator both to the Navier-Stokes equations and to the 
equation for the potential temperature, and exploiting the decomposition 
(|TT| ) (and the analogous for the velocity field) in the advection terms one 
obtains the corresponding filtered equations. For the sake of brevity, we 
report the sole filtered equation for the potential temperature: 

d t 9 = -v a dj - d Q rW (12) 

where Ta are the subgrid turbulence fluxes of virtual temperature (in short 
SGS fluxes). They are related to the resolved-scale field as 

rf = -K H dJ (13) 

K H being the SGS eddy coefficient for heat. A similar expression holds 
for the subgrid turbulence fluxes of momentum (see Moeng, 1984) that are 
defined in terms of the SGS eddy coefficient for momentum (Km)- 

The above two eddy coefficients are related to the velocity scale e' , e' being 
the SGS turbulence energy the equation of which is solved in this model, and 
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to the length scale I = (AxAyAz) 1 ^ 3 (valid for the convective cases) Ax, Ay, 
and Az being the grid mesh spacing in x, y and z. Namely: 



K M = 0.1 I e> 



77V2 



(14) 



K H = 3K M . 



(15) 



4 The simulated convective experiment 

In the present first study, our attention has been focused on the Simulation 
B (hereafter referred to as Sim B) by Moeng and Sullivan (1994). Sim B is a 
buoyancy-dominated flow with a relatively small shear effect, where vigorous 
thermals (see again Fig. 1) set up due to buoyancy force. 
The sole difference of our simulation with respect to the Moeng and Sullivan's 
simulation is the increased spatial resolution, here of 128 3 grid points. 
A preliminary sensitivity test at the lower resolution 96 3 (as in Moeng and 
Sullivan (1994)) did not show significant differences in the results we are 
going to present. Sensitivity tests at higher resolutions are still in progress 
and seem to confirm our preceding conclusion. 

Our choice for a convective boundary layer lies on the fact that, in such 
regimes, dependence of resolved fields on SGS parameterization should be 
very weak, and thus LES strategy appears completely justified. Indeed, in 
convective regimes, SGS motion acts as net energy sinks that drain energy 
from the resolved motion. This is another way to say that energy blows from 
large scales of motion toward the smallest scales and the cumulative (statisti- 
cal) effect of the latter scales can be successfully captured by means of simple 
eddy-diffusivity /viscosity SGS models. Uncertainties eventually present at 
the smallest scales directly affected by SGS parameterizations (that are not 
the concern of the analysis we are going to show) do not propagate upward 
but are promptly diffuse (and thus dissipated) owing to the action of the 
aforementioned eddy-diffusivity/viscosity character of SGS motion. Genuine 
inertial range dynamics can thus develop and, as we shall see, the typical 
features characterizing an inertial range of scales (e.g., rescaling properties 
of statistical objects) to appear. 

The following parameters characterize the Sim B. Geostrophic wind, U g = 
10 m/s; friction velocity, u* = 0.56 m/s; convective velocity, u>* = 2.02 m/s; 
PBL height, Zi = 1030 m; large-eddy turnover time, r* = 510 s; stability 
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parameter, Zi/L = —18 (L being the Monin-Obukov length); potential tem- 
perature flux at the surface, Q* = 0.24 mK/s. 

Moreover, the numerical domain size in the x, y and z directions are L x = 
L y = 5 km and L z = 2 km, respectively; the time step for the numerical 
integration is about 1 s. For details on the simulated experiment, readers 
can refer to Moeng and Sullivan (1994). 

To perform our statistical analysis, we first reached the quasi-steady state. 
It took, as in Moeng and Sullivan (1994), about six large-eddy turnover 
times, r*. After that time, a new simulation has been made for about 37r* 
and the simulated potential temperature field saved at 0.5r* intervals for the 
analysis. Our data set was thus formed by 74 (almost independent) potential 
temperature snapshots. 

Each simulation hour required about 24 computer hours on an Alpha- 
XP1000 workstation. 

5 Results and discussions 

5.1 Statistics of large temperature fluctuations 

Let us start our statistical analysis from the large temperature fluctuations. 
These are controlled by the pdf tails of temperature differences, A r . z 9, and, 
as we are going to show, they are compatible with the laws (§) and (§), that 
means intermittency saturation. 

To show that, it is enough to see whether or not there exist a positive number, 
, such that the quantities P(A r . z 9)r~^°° collapse on the same curve, 
Q, for different values of the separation r. Indeed, as showed in Sec. 0, in 
the presence of saturation the function Q, appearing in @, does not depend 
on r. 

The validity of (§) can be seen in Fig. 2, where the behavior of P(A r . z 8) for 
z/zi = 0.3 and two values of r are shown, z and Zi being the elevation above 
the bottom boundary and the mixed layer height, respectively. In the graph 
(a), P(A r;z 9) is reported without any r-dependent rescaling; in (b) we show 

a^P(A r . z 6)r^^°° for ~ 0.6. The data collapse occurring on the tails of 
the curves of graph (b) is the footprint of intermittency saturation. 
In Figs. 3 and 4 we show the analogous of Fig. 2 but for z/zi = 0.45 and 
z/zi = 0.6. Also in these cases, the exponent giving the data collapse is 

(z) 

Coo ~ 0.6. Similar behaviors have been observed for all z s within the mixed 
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layer. 

As a conclusion, from the evidences of Figs. 2, 3 and 4, it turns out that the 
saturation exponent does not depend on z within the mixed layer. It 
is thus a property of the entire mixed layer and, for this reason, it will be 
simply denoted by Coo- 

Let us now corroborate the scenario of intermittency saturation by look- 
ing at the cumulated probability (|7|). 

For the saturation to occur, such probability has to behave as a power law 
with exponent Coo (see (§)). Such behavior is indeed observed and showed in 
Fig. 4(a) (for z/zi = 0.3 and A = 5 and 5.5) and in Fig. 4(b) (for z/zi = 0.6 
and A = 5 and 5.5). The continuous lines have the slope Coo ~ 0.6 as mea- 
sured from Figs. 2 and 3. The fact that there exist a region of scales, r, where 
that slope is parallel to the slope of the cumulated probabilities means, again, 
intermittency saturation with a unique (i.e. characterizing the whole mixed 
layer) exponent. 

It is worth noting that figures similar to Figs. 5(a) and 5(b) have been ob- 
tained also for smaller values of A, e.g., A = 2.5 and A = 3.5. Population of 
strong events being decreasing as A increases, the above independence on A 
points for the robustness of our statistics. 

In Fig. 6(a) we report, for zjz{ = 0.45, the behaviors of the sixth and 
eight-order structure functions of temperature differences vs the separation 
between points (squares). Stright lines have the slope Coo = 0.6. This is 
a further, direct, evidence of intermittency saturation. To investigate the 
statistical convergence of our sixth and eight-order moments, we reported in 
Fig. 6(b) the bulk contribution to such moments: (A r 8) p P(A r 6), with p = 6 
and p = 8, r/L ~ 7 x 10~ 2 in the inertial range of scales. For comparison 
P(A r 8) is also shown. Note that the maximum contribution to the moments 
six and eight comes from fluctuations A r 8 in the region where P(A r 8) (i.e., 
p = 0) is well resolved, i.e., our statistics appears reliable up to the order 
eight. 

Once Coo is known, we evaluated from (|6]) the unknown function Q. Such 
function is shown in Fig. 7 for two different values of z within the mixed layer. 
Differences among the two curves are evident, signaling that Q contains a 
dependence on the elevation z. Such dependence can be associated to the 
relatively small shear present in our convective simulation (see Sec. f|) . 
Further simulations spanning intermediate ABL (i.e. where both shear and 
buoyancy are important) have however to be performed in order to confirm 
the above last conclusion. 
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It is worth stressing that the scales r at which we observe scaling behaviors 
are always larger than ~ 8 grid-points (i.e. sufficiently far from the scales 
directly affected by SGS parameterizations). Our attention being focused in 
a region sufficiently far from boundaries, this is another point in favor for 
the possible SGS independence of our results. 



5.1.1 A link between geometry and statistics 

Let us now conclude this section with a geometric point of view for the in- 
termittency saturation (see also Celani et al, 2001a). As we shall see, the 
saturation exponent Coo is related to the fractal dimension, Dp, of the set 
hosting the strong temperature fluctuations. 

To do that, let us schematize in a very rough way our strong (i.e. larger 
than some a^) temperature fluctuations in the form of quasi-discontinuities 
(i.e. step functions). Each quasi- discontinuity will define a point (we are 
performing the analysis on planes at constant z) of given coordinate in our 
two-dimensional plane. The ensemble of all points defines the set, S, hosting 
strong temperature fluctuations. Roughly, S is formed by the intersection 
of our two-dimensional plane with the plume interfaces across which strong 
temperature jumps occur. 

A useful indicator to characterize geometrically our set, S, is the fractal 
dimension, Dp, (see, e.g., Frisch (1995) for a presentation oriented toward 
turbulence problems). We briefly recall the standard way to define Dp. 



• Take boxes of side, r, and cover the whole plane at fixed z. Denote 
with N to t the total number of those boxes; 

• Define the function N(r) as the number of boxes containing at least 
one point of S; 

• For r sufficiently small, one expects power law behavior for N(r) in the 
form: N(r) ~ r~ Di % which defines the fractal dimension of S. 

Given the fractal dimension, Dp, it is now easy to compute the probability, 
Prob[\A r . z 9\ > Aer^], of having strong (i.e. larger than some a^) tempera- 
ture jumps within a certain distance r. Indeed, by definition, we have: 

n , n , n . . ,_■>. favorable cases N(r) r~ Dp 9 n , . 
Pro mr .J\ > W»l S posBib[e cases =-±1 ~ — = r>-°' . (16) 
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From (§) and (|S|) the identification = 2 — D F immediately follows. Notice 
that if one does not restrict the attention on the sole planes at constant z, 
but focuses on the whole three dimensional space, the above relation becomes 
Coo = 3 — D' F where D' F is the fractal dimension of the new set S. 

5.2 Statistics of weak temperature fluctuations 

Let us now pass to investigate the statistics of well-mixed regions of the 
temperature field, corresponding to the inner parts of plumes that are likely 
to be present in our CBL (see again Fig. 1). 

In these regions, fluctuations turn out to be very gentle and, as an immediate 
consequence, statistics is expected to be controlled in terms of single-scale 
fluctuations (see the Introduction). The best candidate to characterize, from 
a statistical point of view, weak fluctuations is thus the rescaling form given 
by (§• 

To investigate whether or not our data are compatible with such rescaling, it 
is enough to verify whether there exist a number, aS z \ (a priori dependent 

( 2 ) 

on the elevation z within the mixed layer) such that, looking at r a P(A r]Z 9) 
vs A r . z 9/r a(z) for different values of r, all curves collapse one on the other for 
each value of z. 

Our data support this behavior for the pdf cores (as expected, the rescaling 
(H) holds solely for weak fluctuations), as it can be observed in Fig. 2(c) (for 
z/zi = 0.3), in Fig. 3(c) (for z/ Zi = 0.45) and Fig. 4(c) (for z/z { = 0.6). In 
all cases, the values of is ~ 0.2, that means that aS^ does not depend on 
z. As the exponent Coo, a = or- z > thus characterizes the entire mixed layer. 

6 Conclusions and discussions 

We have characterized, from a statistical point of view, both large and weak 
temperature fluctuations of a convective boundary layer simulated by a large 
eddy simulation model. 

The main results of our study can be summarized as follows. 

• Large temperature fluctuations, occurring across plume interfaces, turn 
out to be strongly intermittent. This is the cause of the observed break 
down of mean field theories a la Kolmogorov, predicting a linear be- 
havior of the scaling exponents, ( p , of the structure functions of tem- 
perature differences, vs the order p. We found, on the contrary, a 
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pdf rescaling which corresponds to a nonlinear shape of ( p vs p, with 
Cp ~^ Coo = const for p large enough. This behavior is named intermit- 
tency saturation, i.e. the strongest violation of dimensional predictions. 
Hence, the concept of 'typical fluctuation' does not make sense: it is 
necessary to specify which part of the pdf of temperature differences is 
sampled for the definition of 'typical fluctuation'. 

• Weak temperature fluctuations, characterizing the inner plume region 
where temperature is extremely well-mixed, have a self-similar charac- 
ter. This amounts to say that, despite the fact that many scales are 
excited in the well-mixed regions, the concept of 'typical fluctuation' 
here makes sense. In this simple rescaling characterizes the pdf 
core, which corresponds to a linear behavior of the curve ( p vs p for 
small p's. The slope of the straight line ( p vs p is a ~ 0.2. 

• Exponents a and Coo appear to be independent on the elevation within 
the mixed layer. They are thus an intrinsic property of the entire mixed 
layer. 

• Statistics and geometry turn out to be intimately related. A simple 
relationship holds indeed between Coo and the fractal dimension, D F , 
of the set hosting the large temperature fluctuations: Coo = d — D F 
where d is the usual dimension of the physical space. 

As for Coo, Dp appears an intrinsic, i.e. ^-independent, property of the 
entire mixed layer. 

It is worth stressing that the present scenario holds also for idealized, 
two-dimensional, models of scalar turbulence both passive (see Frisch et al 
(1999), Celani et al (2000) and Celani et al (2001a)) and active (see Celani 
et al (2001b)), simulated by means of direct numerical simulations. This 
fact naturally points toward the possible generality of the present statistical 
characterization for the entire class of scalar transport problems. 

Finally, let us discuss a possible application of our results within boundary 
layer physics, and, more specifically, in the LES approach. As well known, 
one of the most challenging problem in the LES strategy is to find a proper 
way to describe the dynamical effect of small-scale unresolved motion on the 
resolved large scale dynamics. Recently, new approaches have emerged as 
alternatives to the eddy viscosity and similarity models (see, e.g., Meneveau 
and Katz, 2000). They construct the small-scale unresolved part of a total 
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field (e.g., the velocity field) by extrapolating properties of the (resolved) 
coarse-grained field. A specific form for the subgrid scale field is thus postu- 
lated exploiting scale-invariance (i.e. inertial range scaling behavior) of the 
coarse-grained field. Note that, standard approaches postulate the form of 
the stress tensor rather than the structure of the unresolved field. 
The mathematical tool which permits to perform such an interpolation is 
known as "fractal interpolation" (see Meneveau and Katz (2000) and refer- 
ences therein) where the free parameter of the theory is the fractal dimension 
of the field. 

Up to now, many efforts have been devoted to exploit such strategy for the 
velocity fields. The latter exhibit however a multifractal structure (roughly 
speaking, an infinite set of fractal dimensions characterizes the whole field) 
and thus the fractal dimension parameter is a sort of 'mean field' description. 
The suggestion arising from our results is that the same technique exploited 
for the velocity field appears successfully applicable in the convective case to 
the temperature field as well. Indeed, our results support a fractal structure 
of the temperature field. The situation seems to be even better than that 
for the velocity field. As we have stressed in the preceding sections, solely 
two exponents characterize most of the statistical properties of temperature 
fluctuations. We thus propose to schematize temperature fluctuations as a 
bi-fractal object (i.e., the simplest multifractal object) described by our two 
exponents a and Coo and to generalize "fractal interpolation" to this case. 
To do that, it seems necessary to investigate how the exponents we found 
for the analyzed experiment change by varying, e.g., the weight of buoyancy 
with respect to shear. We are currently working on this point and, as far 
as we remain on convective experiments, it seems that only small variations 
appear. 

It should also be interesting to investigate what happens to the above sce- 
nario in stable stratified boundary layers. In that cases LES approach appears 
more delicate than in a CBL. Observations in field and/or in a wind tunnel 
should be thus necessary to investigate the problem. 
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List of Figures 



1. A typical snapshot of the potential temperature field 6, in the quasi- 
steady state of a convective boundary layer simulated by a Large Eddy 
Simulation with resolution 128 3 . Above: vertical cross-section restricted 
to the mixed layer; below: horizontal cross-section inside the mixed 
layer. Colors are coded according to the intensity of the field: white 
corresponds to large temperature, black to small ones. Plumes and 
well-mixed regions are clearly detectable. 

2. The pdf's P(A r . z 6), for two values of r inside the inertial range of scales 
(dotted lines: r/L = 0.22; continuous line: r/L = 0.11, L being the 
side of the (squared) simulation domain) and z/zi = 0.3, being the 
elevation of the mixed layer top. (a): pdf's are shown without any 
r-dependent rescaling; (b) the pdf is multiplied by the factor a^T~^°° 
with (oo ~ 0.6: the collapse of the curves indicate the asymptotic 
behavior P(A r - z 6) ~ r^ for large A r . z 8, that means saturation of 
temperature scaling exponents (see, and @, (§) and @); (c) pdf's 
are multiplied by the factor r a while A r . z 6 by r~ a : the collapse of 
pdf cores indicates the validity of (|3|) that is equivalent to the linear 
behavior of low-order temperature scaling exponents (see (ffl)). 

3. As in Fig. 2 but for z/z { = 0.45. 

4. As in Fig. 2 but for z/zi — 0.6. 

5. The cumulated probabilities Prob[\A r . z 9\ > \a^] for two values of A 
are shown for (a): zjz\ = 0.3 and (b): z/zi = 0.6. The slope of these 
curves (continuous line) are compatible with the exponent Coo ~ 0.6. 
The error bar on this slope is of the order of 0.1, evaluated by means 
of the local scaling exponents (on half-decade ratios) as customary in 
turbulence data analysis. 

6. (a) Sixth and eight-order structure functions of temperature differences 
vs the separation between points (squares). Stright lines have the slope 
Coo = 0.6. (b) The bulk contribution to the moments p = 6 and p = 8, 
(A r 6) p P(A r 8), with r/L ~ 1.1 x 10 _1 in the inertial range of scales. 
For comparison P(A r 6) (i.e. for p = 0) is also shown. Note that 
the maximum contribution to the moments six and eight comes from 
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fluctuations A r 8 in the region where P(A r 8) is well resolved. This 
proves the reliability of our statistics to compute moments up to the 
order eight. 

7. The function Q defined in (||) is shown for two different values of z 
within the mixed layer: z/Zi = 0.3 (dotted line) and z/zi = 0.6 (con- 
tinuous line). Differences in the shape of these two curves, reveal that 
Q contains a dependence on the elevation, z, within the mixed layer. 
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